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Abstract. Weinberg's dimensional power counting for nuclei is reprised in 
this primer. The role of QCD and chiral symmetry in constructing effective 
Lagrangians is discussed. The two-scale hypothesis of Manohar and Georgi is 
combined with power counting in amplitudes to shed light on the scales of nu- 
clear kinetic and potential energies, the size of strong-interaction coupling con- 
stants, the size of meson-exchange currents, and the relative sizes of JV-nucleon 
forces. Numerous examples are worked out and compared to conventional nu- 
clear models. 

1 Nuclear Perspectives 

What is a nucleus and what are its constituents? This question has no direct 
answer, or even a unique one. Indeed, the proper answer depends on the context 
of the question. Alternatively, the answer depends on the energy scale at which 
we are probing the nucleus. At extremely high energies such as those of CERN 
or RHIC, we might suppose that the appropriate degrees of freedom (d.o.f.) are 
quarks and gluons. At very low energies (such as nuclear ground states), the 
more natural degrees of freedom would be nucleons and possibly pions. After 
all, collisions of a nucleus with low-energy projectiles (~ tens of MeV) eject 
nucleons. This traditional description of a nucleus as interacting nucleons would 
seem at first glance to have little direct connection with QCD, the fundamental 
theory of the strong interactions. The former description in terms of quarks 
and gluons obviously has such a connection, since the beautiful economy and 
symmetry of the theory is manifest in terms of these d.o.fQ, |j. 

A close connection between the two descriptions must exist, but it is murky. 
The traditional picture of nuclei has been very successful, provided that the 
complexities of the short-range nuclear force are resolved by appeal to experi- 
ment. While this works for the nucleon-nucleon (NN) forceQ, we have as yet 
been nearly powerless to resolve the three-nucleon force in this way, although 
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substantial progress is being madcQ. The few-nucleon systems are arguably 
the area of nuclear physics where the most progress has been made on seminal 
problems during the previous decade or sojj|. We expect (and will see) that 
these systems do provide us with substantial information on dynamics, assum- 
ing that we can interpret it. That is the purview of this article. Conversely, 
the scheme that is outlined below has the potential to augment greatly our 
understanding of nuclear dynamics, and particularly the role and appreciation 
of few-nucleon systems (our testing ground) in this endeavor. 

Traditional nuclear physics is the domain of large coupling constants (the 
interaction is strong), which presupposes that perturbation theory cannot con- 
verge. Although we will see below that this is true in part, it is not true in toto, 
and this caveat makes progress possible. Connections to QCD (at least the 
qualitative ones) are easy to make and most are long knownjl], |2j, and follow 
from the chiral symmetry (CS) manifested by QCD. 

• The pion has a very small mass, . 

• The pion is a pseudoscalar particle {J* — CP). 

• The pion is an isovector particle. 

• Chiral symmetry forbids large ttN interactions. 

• There is a large-mass scale, A ~ 1 GeV, associated with QCD. 

The first four items involve the pion and are extremely important, but have 
been known for decades ||. The last item is new|l|, |t], ||, and its consequences 
for nuclear physics are more subtle and will shape the contents of this article. 
In the past, the lack of any apparent dynamical scales in nuclear physics made 
progress in interpreting results rather difficult. 

The (broken) chiral symmetry referred to above is a consequence of (nearly) 
massless quarks in the QCD Lagrangian. This symmetry is reflected in the 
Goldstone mode by (nearly) massless pionsQ. An analogous chiral symmetry 
exists in QED if the electron mass is set to zero; this symmetry produces a 
vanishing amplitude for (relativistic) electrons scattering backwards from the 
nuclear charge distribution. The very small pion mass plays a significant role 
in nuclear physics, producing the longest-range part of the strong force. The 
pseudoscalar nature of the pion guarantees a spin-dependent interaction with 
a nucleon, and this leads to the tensor force, the dominant component^ of 
binding in few-nucleon systems ({V W )/{V) ~ 75%) and possibly in all nuclei. 
The existence of charged pions guarantees them a large role in electromagnetic 
(EM) meson-exchange currents JhJ, which we will treat later. The fourth item 
on the list is critical for a tractable theory of nuclei. It guarantees that a class of 
many-body forces is weak[[ll], |l2| (we will see later that they all are weak). Item 
five allows us to combine phenomenology with principle, and is the organizing 
element of Chiral Perturbation Theory (xPT)@, [| ||, ||, ||], which we will 
treat qualitatively. It is not unreasonable to expect that this subject will become 
increasingly important to nuclear physics in the next decade, and this primer 
was motivated by that supposition. 

As evidence for these views, Table 1.1 shows the (current) results of one 
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of the most important (set of) calculations undertaken in nuclear physics JL6| . 
Using the best available NN force, a weak three- nucleon force (3NF), and no 
four-nucleon force (4NF), the light nuclei (A < 6) were calculated with an 
uncertainty ^ 1%. The agreement with experiment must be considered very 
good. Our goal will be to understand these results in simple terms. 



Table 1.1. Calculated and experimental ground-state energies (in MeV) of few- 
nucleon systems, together with (approximate) dates when they were first accurately 
solved for "realistic" potentials. 



A X(J T ) 


2 H(1+) 


3 H(r) 


4 He(0+) 


5 Hc(§-) 


5 He(i-) 


6 Li(l+) 


Solved 


-1950 


1984 


1987 


1990 


1990 


1995 


Expt. 


-2.22 


-8.48 


-28.3 


-27.2 


-25.8 


-32.0 


Theory 


-2.22 


-8.47(2) 


-28.3(1) 


-26.5(2) 


-25.7(2) 


-32.4(9) 



We summarize this discussion by stating that 

• Chiral Symmetry (manifested in QCD) provides an organizing principle for 
nuclear physics. 



2 Motivation 

In 1979, Weinberg[H suggested that a convenient and simple way to reproduce 
the results of current algebra (and go beyond) was to use a nonrenormalizable, 
phenomenological Lagrangian that manifests chiral symmetry. Such a scheme 
would produce amplitudes of the form: T — E", where E is the energy. This 
was obtained using dimensional analysis and specified v to be an integer deter- 
mined by the type of process and by chiral symmetry. The constraints of CS 
mandate that more complicated mechanisms necessarily have larger values of 
v. Thus, provided that E is smaller than some intrinsic energy scale, A, one 
has a decreasing series (i.e., it is a series in E/A) that is calculable. This se- 
ries is organized in terms of v rather than around unspecified (large) coupling 
constants. For this to work effectively A must be sufficiently large, and more 
complicated mechanisms must not produce smaller values of v. That is, ampli- 
tudes from loops and other products of higher-order perturbation theory (PT) 
should not be larger than lower orders. As noted by Weinberg, the derivative 
structure of his Lagrangian (which enforced CS) guarantees this. Concomitant 
with the nonrenormalizability was the appearance of more and more unknown 
constants at higher orders in PT, which must be determined from experiment 
or calculated from QCD. This feature clearly is a drawback, and one hopes 
that meaningful results can be obtained before the number of terms becomes 
too large for tractability. 

This seminal, but skeletal, argument was developed into a highly success- 
ful program by the Bern group (and others) who introduced fermions. A 
very important addition to the counting argument was made by Manohar 
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and Georgipl, which we will discuss later. Much later (in a series of papers) 
Weinberg[|ll| applied the procedure to nuclei. Different terse versions of the 
counting were constructed in the various papers, while adapting the argument 
to the nuclear environment, a nontrivial achievement. This was extended in 
several ways in the thesis of van Kolck|^|, which is an excellent introduction 
to those aspects of the problem that we will ignore. 

Our task is to reprise the derivation of Weinberg, add expository material 
familiar to nuclear physicists in order to make the results concrete, and finally 
to extract a number of crisp conclusions with import for nuclear physics and 
especially for the few-nucleon field. It can be fairly stated that there is nothing 
in this primer that is not stated or implied in the work of Weinberg, van Kolck, 
and others. Nevertheless, the importance of these ideas and their unfamiliarity 
to our field make an exposition desirable. We will see that many qualitative 
results follow from simple arguments, once the framework has been constructed. 

In order to satisfy the whims of the author, an older alternative form of 
power counting (the "rules of scale" jl8| p^|) will also be separately presented. 
Although this has been in use for decades |20|, it is inappropriate for such com- 
plex mechanisms as loops and, moreover, it has no grounding in chiral symme- 
try. It is, however, rather simple minded, easy to use, and does generate some 
insight into the nuclear physics that is less visible in the more formal approach. 

Finally, we will fashion all of our arguments around Feynman diagrams, so 
we assume that the reader has some familiarity with their struct urcpl[|. Specific 
examples with regularization and renormalization are relegated to Appendices 
B and C. The first example illustrates a number of currently popular techniques 
and approaches, while the second treats a simple two-nucleon problem. An 
exposition of nuclear matrix elements in momentum space is extremely helpful 
in interpreting the power counting in nuclei, and is relegated to Appendix A. 
Otherwise, we will attempt to motivate the physics as we proceed. 

We will show the following results at various points in the primer: 

★ The kinetic energy /nucleon scales as Q 2 /A, where Q is an effective mo- 
mentum in the nucleus. 

★ The "intrinsic" potential energy/pair scales as -Q 3 / f^A, where / w is the 
pion-decay constant. 

★ The cancellation of these two comparable energies allows nuclei to be 
weakly-bound systems ((T) ~ — (V) and binding energy <C mass). 

★ Large strong-interaction coupling constants are caused by the mismatch 
between the /„■ and A scales. 

★ A^-nucleon forces decrease in strength as N increases. This powerful result 
authenticates decades of nuclear physics supposition and phenomenology. 

★ Increasingly more complicated forces contribute more weakly. 

★ Arbitrary processes contribute to the energy as Q" ', where v = 1 + 2(n c + 
L) + A, with L the number of loops, n c a topological parameter, and A > 
reflects the complexity of the interaction. This is Weinberg power counting. 



3 Rules of Scale - Nuclear Methodology 
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A type of primitive power counting has been in use for two decades@ ^ |o). 
The "rules of scale" were developed as a way to control expansions that arise 
in nonrelativistic treatments of meson-exchange currents. Specifically, expan- 
sions in powers of 1 /c 2 were organized by noting that the velocity of a typical 
nucleon (with mass Mjq and momentum p) is v = p/Mjsr, and consequently 
v/c ~ p/Mnc. Thus, counting powers of 1/M/v in an expansion is equivalent 
to counting powers of 1/c. It was noted that nuclei are weakly bound, implying 
that the kinetic and potential energies satisfy (T) ~ — (V). Hence, V should be 
counted as 1/Mjv, since V ~ T ~ p 2 /2Mjv- This scheme works well at "tree" 
level (i.e., those mechanisms that don't involve loops, which were typically of 
interest in those times) and, more importantly, it works for the current continu- 
ity equation, which is essential for EM processes. That is, consistency in that 
equation demands that energies of all types be treated on the same footing. 
Conversely, the scheme has no predictive power for loops. 

In order to progress further, we need some indication that v/c ^ 1. Oth- 
erwise, the Schrodinger approach is meaningless. We can use the uncertainty 
principle and the radii [^3] of the few- nucleon systems (~ 1.5 - 2.0 fm) to pro- 
duce: pc ~ fic/R ~ 100-150 MeV ~ m,c 2 , where the last relation is a conve- 
nient mnemonic, rather than a statement of principle, since m w (unlike f w , for 
example) vanishes in the chiral limit. This leads to 

v/c ~ p/M N c ~ m^/Mjv - 10 - 20% , (3.1) 

and (v/c) 2 is satisfyingly small on average. In order to conform to a more 
standard notation, we will use Q (~ m OT ) rather than p. We will see later that 
in processes involving pions this is indeed a typical scale 

At this point a little history is instructive. There are at least three ways of 
organizing the interactions of pions with nucleons. This familiar argument has 
been distracting and ultimately unproductive (the author has been guilty of 
thisp9[ |20| , |25|), but it illustrates a number of useful points. Parity and time- 
reversal-invariance arguments allow the Lagrangian for a single pion-nucleon 
interaction to be either 

C wNN = -iGN~/ 5 T-irN, (3.2) 

or 

£' wNN = -J-N l5 p(T-?r)N, (3.3) 

where / = Gm^ /2M^ . These two forms are usually referred to as PS (pseu- 
doscalar) and PV (pseudovector) . An alternative form is introduced by using 
the Goldberger-Treiman ]26| relation: 

where G ~ 13, qa — 1-26 is the axial-vector coupling constant and f n ~ 92.4 
MeV is the pion-decay constant. This remarkable relationship (the author's 
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favorite in strong-interaction physics) relates strong interactions on the left 
side to weak interactions on the right, and is violated at the level of 2% [^?J by 
chiral-symmetry breaking (the left side is larger) . Thus we can also rewrite Eq. 
(3.3) by replacing f /m^ with 17^4/2/^, generating the third and preferred form 

C , ; NN ^-^N 1 rJ(t.7v)N, (3.5) 

where we have written t = r/2 (a very common practice). 

Note the very different sizes of the interactions in Eqs. (3.2) and (3.3): 
G ~ 13, while / ~ 1. These two forms are equivalent on-shell, but differ dra- 
matically off-shell. We will see later that the off-shell behavior is controlled by 
chiral symmetry through terms we have not explicitly written. In particular 
the large iViV- "pair" terms implicit in the first form are exactly cancelled, a 
phenomenological result known historically as "pair suppression" . These pair 
terms have a very large dimensionless strength, G 2 m 7V /Mpf ~ 25, producing 
many-body forces that would make the tractability of nuclear physics prob- 
lematic. The large size of G is typical of the strong interactions, while / is 
much more modest. Even more modest is the effective coupling constant for 
the one-pion-exchange potential (OPEP) shown in Fig. (la). This amplitude, 
when converted to configuration space, generates a factor of 1 /Air, whose precise 
form depends on the fact that we live in three space dimensions (see Appendix 
B): 

f d 3 q e z Q ' r ma 
J (2?r) 3 q 2 + ml ~ Air 



(3.6) 




(a) (b) (c) (d) 

Figure 1. One- and two-pion-exchange contributions to the nuclear force. Solid lines 
are nucleons, while dashed lines are pions. 



Rendering OPEPH into a dimensionless radial function of x — m T r and 

3 2 

an overall dimensionful factor, we have V„(r) = ^Jj#-V(x). Obviously some 
care is required in this exercise, since Air is dimensionless. It is, however, a large 
number and we extract it. We expect (V(x)) ~ 1 since (1) = 1, where 1 is the 
unit operator (this is our normalization convention) and because correlations 
with a length ~ l/m^ are introduced into the wave function by V(x). This 
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implies that values of x near 1 are very important (which can be verified in 
Fig. (7)) and our functions near x = 1 are 0(1) [cf., Eq. (3.6)]. The net result 
is that we expect a pair of nucleons exchanging a pion to generate an energy 
~ m^g\i A-n f% ~ 20 MeV. Anticipating a result from Section 5 and Appendix 
B, we note that Anf^ ~ 1200 MeV ~ A, where A was introduced in Section 
1. Recalling Q ~ m^, we obtain (OPEP) ~ -^j-Cy per pair, where Cy is 
the product of all the dimensionless factors ~ 1 (by supposition). We have 
already argued that (T) ~ — ^ T per nucleon, where Ct ~ 1 and we have used 

Mjy ~ A. For the triton, 15 MeV/pair and (T) - 15 McV/nucleon[|4). 

For the a particle both energies are somewhat larger. In heavier systems, where 
the Pauli Principle begins to play a large role, these numbers (particularly the 
potential energy) cannot be correct. Nevertheless, the "intrinsic" scale for these 
quantities is 

(T) /nucleon- Q 2 /A, (3.7) 

' (3 ' 8) 

confirming that kinetic and potential energies (and quantities derived from 
them, such as impulse-approximation and meson-exchange EM currents) should 
be counted similarly (numerically, Q ~ /„■), in the absence of large dimension- 
less coefficients that would skew the similarity. Thus, for OPEP and light nuclei 
we see a confirmation of the weak-binding hypothesis that formed the basis for 
the rules of scale. Moreover, the numerical results based on calculations with 
realistic potentials are consistent with these scales. 

Two other potentials have been calculated that can be similarly analyzed. 
Figures (lb) and (lc) show various parts of the two-pion-exchange potential 
(TPEP). These Fcynman graphs can be separated into two components: an 
intrinsic TPEP and the simple iteration of OPEP (which happens automatically 
whenever the Schrodingcr equation is solved) . The "intrinsic" size of the former 

Two-pion-exchange three-nucleon forces (2ttE3NF) can be similarly 
analyzedjjD (see Fig. (Id)): 

m 6 O 6 3 

{2nEiNF) ~ Msfrurz ~ ^ ~ i MeV/triplet ■ (3 - 10) 

In counting arguments of this type, we have avoided factors of ^u, 2, 1/2, etc., 
and hope that they average out ~ 1. In most cases, they do (within an overall 
factor of 2). In Sections 6, 8, and 9, we will derive a very simple formula 
(Weinberg's formula) that reproduces all of these results by power counting (in 
Q and I /A) with almost no effort. 

The only remaining problem is Fig. (lb). This diagram, unlike a typical 
case, does not have momenta — Q ~ flowing through every propagator. 
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Indeed, after the first interaction the propagator is typically (E — T) _1 ~ 
(p 2 /2Mat) _1 ~ Mff/Q 2 . Graphs of this type are called "reducible". The small- 
er or "infrared" singularity enhances the graph by a factor ~ M/v/m^ com- 
pared to a typical case. This singularity is also the origin of the "ambiguity 
problem" pi ^0|, , which reflects the fact that a potential is an off-shell am- 
plitude, and hence is not unique. The dichotomy between these two types of 
graphs (reducible and irreducible) greatly reinforces the case for the following 
calculational scenario: 

o Calculate irreducible graphs (where simple power counting works) and 
define this as the nuclear potential. 

o Solve the Schrodinger equation with that potential - the infrared enhance- 
ments will happen automatically. 

On the basis of several examples (this does not in any sense constitute a 
proof, which will follow in Section 6) we find that 

• Various nuclear energies behave as Q v / A v ~ n "~ l for some values of v 
and n Cl with more complicated mechanisms generating larger values of v. We 
have used A ~ 47r/ ff ~ Mn in this simple counting exercise. 



4 Effective Interactions 



The basic idea is a simple one. Strong-interaction physics can be divided into 
two parts: short-range (high-energy) parts and long-range (low-energy) parts. 
We hope that the former is not dominant, but it could be. This is the most 
difficult domain of nuclear physics, where phenomenology reigns supreme. The 
long-range physics is of two (nonexclusive) types: (1) from pion-range physics, 
and (2) from iterations of the nuclear potential, as we saw in the previous 
section. Even though convergence requirements restrict us to the low-energy 
regime, it is essential that the pion degrees of freedom be treated explicitly; 
they have a strong energy dependence even at low energy. We also saw that 
the infrared singularity that results from iterating the potential plays a very 
important role in nuclear physics. 
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Figure 2. Unflavored meson spectrum [plus Zi-isobar excitation energy] below 1 GeV. 



The short-range regime is where all the complexity of QCD resides. Con- 
sider, for example, the unflavored meson spectrum below 1 GeV, plotted in 
Fig. (2). Only two mesons lie below 770 MeV: it and rj. In an SU(3) (rather 
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than SU(2)) treatment of the chiral symmetry (i.e., strangeness is included), we 
would treat the (tt, 77, K )-meson octet together. Although the f7-meson contri- 
bution to the nuclear force is not particularly important, in an SU(2) approach 
its low mass would be problematic for the formalism. The other (heavy) mesons 
have masses > 770 MeV. We can therefore imagine a scheme where the effect 
in appropriate spin and isospin channels of all single and multiple heavy-meson 
exchanges, loops, ... is "frozen out"[|l5). That is, the effects of all possible 
short-range mechanisms are lumped together without regard to their origin. 
This is accomplished by noting that (for low energies) a derivative expansion 
about the zero-range limit is formally possible. Thus A ~ m p ~ 1 GeV is also 
the boundary between the short-range and long-range physics. 



(0 



X 




(a) 



(b) 



(c) 



(d) 



(e) 



Figure 3. ui- meson contributions in (a-c) to short-range NN force in (d), plus an 
additional 3NF in (e). Solid lines are nucleons, while wavy lines are ui mesons. 



Figure (3) illustrates this approach, where a few selected mechanisms in- 
volving lu mesons are shown. All are subsumed in (3d), which is the sum of 
various zero-range interactions that can include derivatives. As an example, we 
treat Fig. (3a), for which the S-matrix is given by 



5' 



(-ui7 M ui)(u 2 7' t '"2) 



(4.1) 



where we follow the conventions and metric of Ref. We have dropped 
normalization factors and the overall four-momentum-conserving factor, 
(2n) 4 5 4 (Pi + P2 — P[ — P%), for this two-nucleon cluster. The strong coupling 
of the lu (with mass m w ) to the nucleon is reflected by the coupling constant, 
g u , the nucleons (with spinors, u) have been arbitrarily labeled "1" and "2", 
and g M is the momentum transferred in the collision {q 2 < 0). We now expand 
the denominator for — q 2 <C m 2 . 



1 



1 



(4.2) 



and define each term in the series as a new effective interaction (Lagrangian) 

S ~ i (H) tei 7/t Ui )(U 2 ^U 2 )(1 + At + • ••) , (4.3) 
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leading to (see Ref. fl22[) 

£ e //=4//+4//+-- - . ( 4 - 4 ) 

We have defined the effective Lagrangian so that the usual Feynman rules 
produce, term-by-term, the same S-matrix. We have added a superscript to 
C, U A" = 0, 2,4, . . ., that serves the same purpose as Weinberg's 'V that we 
introduced earlier. It classifies the interaction here according to powers of q: 
we will develop the rules later for this notation. More complicated short-range 
forces are possible, as indicated by the three-nucleon force in Fig. (3e). 

We have used the fact that the u> meson is sufficiently heavy that (for small 
q 2 ) it propagates only a very short distance between emission and reabsorption 
(a la the uncertainty principle) [p8| . Since we have assumed a particular mecha- 
nism, we have developed a model for the strength of the zero-range interaction. 
Unfortunately, a vast array of mechanisms including but not restricted to Figs. 
(3a) - (3c) (viz., loops to arbitrary order, physics not involving the to . . .) also 
contribute. Thus, this coefficient (and others with A > 0) is unknown and must 
be determined phenomenologically or calculated directly from QCD. 

Two more aspects are important and somewhat controversial. In addition to 
mesons being "exchanged" , we also have to consider nucleon resonances. Figure 
(2) shows the Z\-isobar excitation energy superimposed on the meson spectrum. 
The low excitation energy of the A has long stimulated the imagination of nu- 
clear theorists, and mechanisms proposed for almost every nuclear effect will 
typically include a A somewhere. The uncertainty principle states that the 
virtual excitation of a A scales inversely with the A energy (i.e., as 1/Ea)- Re- 
cent work suggests that in some regimes (viz., three-nucleon forces and meson- 
exchange currents (29|) Ea behaves more as A ~ 1 GeV than Ma — Mn ~ 300 
MeV. In other regimes (threshold pp — > ppir , for example |3fj|]) the opposite is 
more likely. In constructing effective interactions, some theorists [[l2) routinely 
include the A in the active Hilbert space, while others do not. Care should 
be taken when assuming that the A is not very important. We will, however, 
ignore the A in what follows solely for reasons of simplicity. 

The second aspect concerns relativity. The simple example involving the 
to meson respected special relativity, including the derivative terms (~ q 2 ). 
More complex mechanisms would lead to terms proportional to the overall 
momentum p M of a nucleon, rather than the momentum transfer, q^. While 
the latter is typically not large, p M = (E,p), and E ~ Mjv ~ A is large. This 
would lead to a series of terms that would all be the same size (recall that we 
have a series in E/A). The problem can be formally eliminated by performing 
a 1/Mjv expansion (e.g., a Foldy-Wouthuysen reduction) 20 ^fj), where the 



effect of NN "pair" terms is projected out, and one treats only the positive- 
energy spectrum. This "freezing out" of degrees of freedom is familiar to nuclear 
physicists in the Feshbach (P-Q) reaction theory[^2[, where the Hilbert space is 
compacted from the normal to a restricted size, leading to much more complex 



11 



operators (which must reproduce the physics of the larger Hilbert space in the 
smaller one). Covariance is lost, but one retains an expansion in Q, rather than 
introducing terms ~ Mn at every order. This alternative formalism is given 
the generic name "heavy-baryon" and is currently the method of choice p^, 
. This presupposes that a nonrelativistic treatment of the nuclear physics 
is appropriate, which clearly holds for the few-nucleon systems. We give an 
example of a nonrelativistic loop calculation in Appendix B, and compare it 
to a relativistic one. Appendix C treats a simple nonrelativistic two-nucleon 
model. 

We summarize this section by noting that: 

• The complexities of strong-interaction physics can be divided into long- 
range and short-range parts, with pion d.o.f. dominating the former and the 
latter condensed into zero-range interactions of unknown size. 

• Effective interactions are constructed by "freezing out" the short-range 
degrees of freedom (mesons, resonances, NN pairs, . . .), leading to structurally 
more complex interactions that (hopefully) are easier to treat at low energies. 

• The A isobar can be included in the active Hilbert space (with the pions) 
or not, depending on the problem to be solved. 

• A nonrelativistic, nonrenormalizable field theory is possible to construct 
and to use. 



5 Dimensional Power Counting in Lagrangians 

We begin by considering the dimensions of various quantities in order to assess 
the scales of the strong interactions. Our approach in this section will be to 
motivate rather than to attempt rigor. The interested reader should consult 
Refs. 0, H for a more sophisticated approach. We will also formulate the 
arguments supposing a relativistic field theory. Following the discussion of the 
previous section, they also apply mutatis mutandis to a nonrelativistic "heavy- 
baryon" field theory. 

We denote a power v of energy by \E V ] (equivalent to length~ v ), and recall 
that a Lagrangian has dimension 4 (i.e., behaves as [E A ]). The Lagrangian for 
a system of pions and nucleons interacting through derivatives (d^) acting on 
any field can be written in the schematic form 

C = C free + AC, (5.1) 
C free = N(i 7lM d» ~ M N )N + i(^7r) 2 - im^Tr 2 , (5.2) 



\ b J \ c ' \ d 

The combination NN has dimension [E 3 ] , while it and d* 1 (as well as ra, and 
Mm) have dimension [E], and AC of course has dimension [E 4 ]. 

Our first assumption is that only the scales and A occur in a,b,c,d 
and also that AC should reproduce Cf ree . Clearly, c must be f n , since that 
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quantity sets the scale for the pion field (cf., PCAC§ § § |l|). Settin S ^ = 
and m = n = 2 gives a/c 2 d 2 ~ 1, while m = 0, = 1, and n = 0, 1 give a/6 ~ A 
and a/bd ~ 1 (where we have substituted Mjv ~ A). These three equations can 
be uniquely solved to produce d ~ A, a ~ f^A 2 , and 6 ~ /^yl and thus 

^ £ = U J - M ' (5 - 4) 

where cg mn is dimensionless, and we have used the form of Cfree to incorporate 
the chiral-symmetry-breaking pion-mass term into Eq. (5.4): a pion mass (either 
explicit or implicit) is treated as a derivative. All Dirac matrices, nucleon isospin 
operators, etc. have been ignored and are indicated by the dots. 

This lovely formula has profound implications if we invoke chiral 
symmetry^], [35]. No unique representation exists for incorporating chiral sym- 
metry (cf., PS vs. PV forms), but a representation always exists in terms of 
covariant derivatives jj], |[ |[ [l2), which has only increasing powers of (1 /A). 
According to Eq. (5.4), these powers have the form (1/A) A , where 

A = n + £-2. (5.5) 

The minimum value for pions is from Eq. (5.2) and for nucleons is —1 for 
each of the two (largely) cancelling terms in Eq. (5.2) for free nucleons (i.e., 
joE— Mm)- This cancellation highlights the problem with derivatives of nucleon 
fields that we raised earlier: a slowly-moving nucleon actually has a kinetic 

energy, Cf ree = —N^^jj^N, or A = 1. The PV (derivative-representation) 
Lagrangian has A = 0. These examples motivate the rigorous result [|| that 
the derivative-representation chiral constraint has a form guaranteeing only 
positive powers of 1/A in Eq. (5.4): 

A>0. (5.6) 

The PS form in Eq. (3.2) has A = — 1, which tells us that additional terms 
in the Lagrangian are required in order to construct a proper representation 
of chiral symmetry. Indeed, it was the neglect of the latter terms that led to 
the ad hoc "pair suppression" mechanism many years ago. This mechanism 
is automatic for chiral models or theories, as noted in many places. The use 
of a chiral representation corresponding to Eq. (3.2) requires an additional 
(leading-order) "pair" -killing interaction of the formp9|: 

AC ps = T^-Ntt 2 N . (5.7) 

We strongly argue that derivative forms be used, because the chiral constraint, 
Eq. (5.6), applies term-by-term. 

The second assumption is that a reasonable theory should havepGj] 

(5.8) 
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or that the theory is "natural". This is also called naive dimensional power 
counting (NDPC)@, ||. Clearly, if the scaling in Eq. (5.4) is a figment of 
our imaginations, the values of C£ mn will jump all over. Verifying naturalness 
validates Eq. (5.4). We also note that if a symmetry exists, ci mn could be 
vanishingly small, but if the scaling hypothesis holds, very small coefficients 
would otherwise be just as unlikely as very large ones. 

The third assumption is that vacuum fluctuations do not alter the scales 
that we introduced. It is shown in Refs. |l], Q that it is necessary for loop 
integrals to be cut off at energies 

A < A-kU (5.9) 

in order for the structure in Eq. (5.4) to be preserved, and we indeed have used 
471"/^ ~ A interchangeably in Section 3 (see also Appendix B). 

At this point, examples will serve us best in assessing the concept of "nat- 
ural", and to illustrate the use of Eq. (5.4). Evaluating that equation for 
£ = m = n = 1(A = 0) and comparing it to Eq. (3.5) produces 

c*~g A ~ 1.26. (5.10) 

Using r/2 instead of t would have produced the equally satisfactory 

The zero-range NN interaction produced by an w-meson exchange was 
derived earlier (Eq. (4.5)) and can be compared to the case I = 2, m = n = 0, 
yielding c u / f% and a natural coupling constant 

~l-75, (5.11) 

using the numerical entries in Table (A. 2) of Ref. p7[ . A famous example is 
the KSFR relation [|38|, which states 



9p_ J_ 



(5.12) 



and corresponds to a zero-range Lagrangian: ^^(Nj^tN) 2 = -Aj^Nj^tN) 2 
or c p ~ 0.25. In general a heavy-meson exchange with mass m x and coupling 
constant (to the nucleon) g x corresponds to a zero-range Lagrangian with a 

2 

coefficient, ~ |f, implying that 

9x~ ^V2^~ T - 10, (5.13) 

J7T Jit 

using m x ~ A, and \J2c x ~ 1. This illustrates why strong- interaction coupling 
constants are large dimensionless numbers (~ 10): the scales A and f„ are very 
different. 

The next example of this type is a caution. What if there were a mechanism 
such as an NN force component of normal size mediated by the exchange of a 
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light meson (m -C A)? This would correspond to a large value of c x , and would 
likely lead in perturbation theory to a growing series whose structure violates 
Eq. (5.4). This is analogous to the intruder-state problem in nuclear-structure 
theory [^9), where a low-energy state badly affects convergence of perturbation 
theory calculations. Another way of saying the same thing is that there would 
be another scale in the problem, which would violate the assumptions that led 
to Eq. (5.4). 

We can also extend the NDPC to other forces. One isospin violation (IV) 
mechanism [[12) is parameterized by the up-down quark-mass difference — 
m u = t{md + m u ), where e = [m,d — m u )/(md + m u ) ~ 0.3. The factor of 
m,d + m u is proportional to to^|| and Eq. (5.4) can still be used, provided that 
we identify [gjj 

cjv ~ec, (5-14) 

where c is 0(1) and n > 2. That is, Eq. (5.4) has the C£ mn replaced by eci mn , 
and there are at least two implicit powers of m,. Similarly, parity- violating 
forces have c replaced bypcf 

— 

where Gf is the Fermi constant and c is 0(1). 

Another application is the nuclear density, which for nuclear matter has 
the "empirical" value Q pnm = 0.153 far 3 ~ 1.5/ 3 . The generic Lagrangian 
contains powers of ( ^^i\ N ^j , the numerator of which is essentially the nuclear 
density (N'N). At nuclear-matter density, the Lagrangian series (in £) is there- 
fore geometric in -jr^f = 1.5/^/ A ~ 1/7, which is comfortingly small, so that 
the expansion presumably works at normal nuclear density. 

Our final example takes us outside the realm of few-nucleon systems. One 
might ask whether there have been any calculations performed using a La- 
grangian of zero-range form. Although this has not been done in few-nucleon 
systems, there does exist one comprehensive Dirac-Hartree calculation jl2| for 
a set of 57 nuclei using such forms (but without pion d.o.f). The various terms 
in the Lagrangian can be represented schematically by a(NN) 2 + (3(NN) 3 + 
j(NN) 4 + 8(d^ NN) 2 , where various isospin and Dirac matrices sit between N 
and N. A total of nine coupling constants were adjusted to fit the energies and 
radii of three representative nuclei. The results are among the best ever ob- 
tained for such a comprehensive set of nuclei. If one uses Eq. (5.4) with m = 
and various values of £ and n, and compares to the expressions in Ref. j|2| |]|, 
one obtains the results in Table (5.1). 

The unsealed coupling constants (whose subscripts refer to Dirac and isospin 
matrices) span 13 orders of magnitude. Most of the scaled c's that result are 
numbers near one. The average of the 7 terms is also natural. The uncomfort- 
ably large difference is possibly due to a lack of sensitivity in that quantity to 
the data. If the tiny value of ars is not an artifact of the fitting process or of 
the neglect of pion-range physics (unknown at this time) it presupposes a sym- 
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Table 5.1. Optimized coupling constants for the Lagrangian of Refs. |42j, [l3[ and 
corresponding dimensionless coefficients and chiral expansion order. 



Constant 


Magnitude 


Dimension 


4mn 


Order 


a s 


-4.508xl(T 4 


MeV~ 2 


-1.93 


A 


OtTS 


7.403xl(T 7 


MeV" 2 


0.013 


A 


ay 


3.427xl0~ 4 


MeV~ 2 


1.47 


A 


axv 


3.257xl(T 5 


MeV~ 2 


0.56 


A 


Ps 


l.HOxlCT 11 


MeV" 5 


0.27 


A- 1 


Is 


5.735xl0~ 17 


MeV" 8 


8.98 


A- 2 


Iv 


-4.389xl0~ 17 


MeV" 8 


-6.87 


A- 2 


Ss 


-4.239xl0~ 10 


MeV" 4 


-1.81 


A- 2 


Sv 


-1.144xl0" 10 


MeV~ 4 


-0.49 


A- 2 



metry of some kind. Nevertheless, scaling as predicted by Eq. (5.4) is obvious. 
Improvements in the quality of the many-body techniques are not expected to 
alter this conclusion 0|, but will alter each number. 
We summarize this section by noting that: 

• A two-scale hypothesis (/„■ and A) for the dimensional factors in the 
effective Lagrangian (plus chiral symmetry) suggests a convergent expansion of 
the Lagrangian series for normal nuclear conditions. 

• This expansion has been organized so that vacuum fluctuations do not 
alter the form. 

• The dimensionless coefficients in the Lagrangian are "natural" if they are 
0(1), which we illustrated by several examples. 

• Chiral symmetry guarantees that the large scale [A ~ 4^/^ ~ m p . . .) 
does not occur with positive powers. 

6 Power Counting in Amplitudes 

Power counting in amplitudes is a straightforward exercise, but somewhat te- 
dious. It also is dependent on the environment and on what one chooses to 
emphasize. The result is worth the effort, however, since most of our previous 
results are subsumed by the very simple final forms. We work with nuclei, bound 
states of A nucleons. Bound states are forever interacting, and one cannot sim- 
ply separate interaction diagrams into connected and disconnected parts and 
discard the latter. Eventually each nucleon interacts and shares energy and 
momentum with the others. In normal (textbook) applications disconnected 
diagrams don't contribute. We emphasize that the sharing of momentum dom- 
inates the systematics in a nucleus. The momentum that one nucleon takes 
from the "bank" is not available to the others. 



16 



The basic idea is to use our knowledge of the dimensionality of propagators, 
phase space factors, vertices, and delta functions to find the dimensionality 
of an amplitude, after eliminating all the coupling constants. Given this and 
knowledge of how the other scales (viz., f v and A) come into the problem, the 
coupling constants can be put back and a complete scheme can be constructed. 
We note that finding a momentum behavior ~ Q v does not necessarily imply a 
final (Q/A) v behavior. We also choose to work in configuration space, and this 
choice means that, while counting powers of momenta, the effect of additional 
phase space factors (besides those in loops) needed to convert to that space 
must also be incorporated, and this has been done in the derivation leading 
to Eq. (6.12). We will give examples later that spell out the differences, and 
Appendix A documents the various factors that normally arise. Finally, we 
postpone dealing with reducible (infrared-singular) diagrams until later. We 
first deal with the amplitudes, and then we will worry about the coupling 
constants. 

Our first concern is what we should calculate and what rules apply. In 
general [pp| , the energy shift, AE, in an interacting system is given in terms of 
the S-matrix (see Appendix B) by S f i = 6 f i — i AE [(2tt) 4 (5 4 (P/ - Pi)}. The 
dfi factor for N noninteracting particles is the product of three- momentum- 
conserving (5-functions of each particle (see Eq. (6.35) of Ref. and behaves 
as [i?~ 3Ar ], while the energy shift (which is also the T-matrix) therefore behaves 
as [E 4 ~ 3N ]. For a cluster of two nucleons (N = 2) this has the dimensions [-E -2 ], 
in agreement with Eq. (4.1) (the spinors and coupling constant in that equation 
are dimcnsionless). That form, as argued above, lacks the phase-space factor 
~ [E 3 ] needed to convert to configuration space. Together these factors produce 
~ [E] , the correct dimension for an energy. We therefore use the same rules for 
an energy shift as for an S-matrix, stripping off the four-momentum-conserving 
(5-function (if there is a single cluster) to produce AE. Note that if there are C 
separate clusters, there will be (C — 1) cluster (5-functions remaining to treat. 

Our second concern is that there are many possible quantities that specify 
how momentum and energy flow in a given diagram. The key is to decide 
which ones to keep and which to eliminate. Typically one eliminates internal 
variables and keeps external ones. This is not enough for specificity and one 
additional choice remains. That choice will be made in such a way that the 
chiral constraint, Eq. (5.6), can be implemented by inspection. We also will 
calculate for fixed A, and eliminate at the end some extraneous factors that 
are dependent on our wave-function normalization scheme and depend only 
on A. Finally, in keeping with common practice, we work in D space-time 
dimensions (D = 4). For simplicity we restrict ourselves to no external bosons 
(interacting with our nucleus) and only calculate the energy shift. 

The first set of variables specifies what happens at individual vertices or 
interaction points, and we assume that there is at least one vertex. Figure (4) 
contains several examples from previous sections plus an additional one. At the 
jth ver t ex j n a Feynman diagram we have 

o bi bosons (i.e., pions in our case) entering or leaving; 
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(a) (b) (c) (d) 

Figure 4. Elements of Lagrangian (i.e., vertices) for a system of pions and nucleons. 
Solid lines are nucleons, while dashed lines are pions. 

o fi fermions (i.e., nucleons) entering or leaving; (6-1) 
o di derivatives acting at the vertex; 
oAi = di + fill - 2 > , 

where the last definition anticipates our final form and corresponds exactly to 
Eq. (5.5), since I — fi/2 and n = di. In Fig. (4a) we have b = 1,/ = 2. For 
the vertex corresponding to Eq. (3.2), d = and A — — 1, while d — 1 and 
A = for Eq. (3.3). Figure (4b) has / = 2,6 = 2, while Eq. (5.7) has d = 
and A = — 1, as we discussed before. Figure (4d) is the Weinberg four-pion 
interaction |8[ [Tl| and has 2 derivatives producing / = 0, 6 = 4, d = 2, and 
A = 0. Recall that the Z\ = — 1 terms separately violate chiral symmetry and 
cancellations must occur between them. Figure (4c) and Eq. (4.5) have / = 4 
and d = b = A = 0. 

The second set of variables will completely specify the process: 

o v = dimensionality of the amplitude without coupling constants {Q u )\ 

o D = number of space-time dimensions; 

o A = number of nucleons (which is conserved); 

o L = number of loops (i.e., internal phase-space integrals) (> 0); (6.2) 
o n c — number of nucleons interacting with at least one other minus the 
number of clusters with at least two interacting nucleons (> 0); 
L = n c + L(> 0) ; 

o A = J2i Ai (> 0) (i.e., sum over all vertices). 

Figure (5) shows a representative case with the requisite complexity. It 
has A = 6 nucleons propagating from some initial time to some final time 
(indicated by cross-hatching for emphasis). The process can be divided into 
two interacting clusters and a single noninteracting nucleon, producing n c — 3. 
There is a closed loop on the left-most nucleon, so that L = 1 and L = 4. 
The very important topological parameter n c specifies the complexity of the 
interaction scenario and will determine the relative importance of iV-nucleon 
forces. 
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Figure 5. System of 6 nucleons interacting via pions. The cross-hatched (for empha- 
sis) regions depict initial and final times. Solid lines are nucleons, while dashed lines 
are pions. 



A third set of variables is required to set up the counting, but will not ap- 
pear in the final result: 

o En = number of interacting nucleons; 
o Xn = number of noninteracting nucleons; 

o C = number of interacting clusters of nucleons; (6-3) 
o Ib = number of boson (pion) propagators (i.e., internal lines); 
o Ip = number of fcrmion (nuclcon) propagators (i.e., internal lines); 
o V = number of vertices (i.e., ^V). 

Note that a cluster must contain at least one vertex, but can have any number 
(including only one) of connected nucleons. Moreover, the nucleon lines touch- 
ing the cross-hatched regions are not propagators. For the case specified by Fig. 
(5) we have C = 2, E N = 5, X N = 1, V = 6, I B = 4, and I F = 1. 

We have several obvious relationships between previously defined variables 



A = X N + E N , (6.4) 
n c = E N -C, (6.5) 

i 



Note that Eq. (6.5) does not appear to correspond precisely to Eq. (6.2). A 
cluster of one nucleon interacting only with itself will not contribute to n c , and 
the two definitions are thus equivalent. Using the fact that two boson (pion) or 
fermion (nucleon) fields are needed to make each boson or fermion propagator, 
we have 

5> = 2/ B , (6.7) 
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fi = Hp + 2E N , (6.8) 

i 

where by definition only interacting nucleons connect to a vertex. The number 
of loops in a diagram can be counted by noting that in constructing a diagram 
there is a four-dimensional phase-space integral associated with each propa- 
gator. Feynman diagrams conserve four-momentum at each vertex, and these 
constraints (V of them) eliminate the integrals. This overcounts, however, since 
each cluster has an overall ^-function (four-momentum) constraint (see the dis- 
cussion below Eq. (4.1)), and there are C of them. Thus, the number of loops 
is 

L = I B + I F -V + C , (6.9) 
and associated with each loop is a four-(i.e., D-)dimensional phase-space inte- 
gral J (2*) l> ~ [E D \- Figure (5) has one loop, by inspection or by using Eq. 
(6.9). Each noninteracting nucleon contributes a three- momentum (i.e., D — 1) 
((-function ~ [E 1 ^ ] (see Appendix A). 

Finally, the momentum or energy dimensionality of any diagram can be 
determined by counting: 
o phase-space factors ~ [E D ] in loops; 
o boson (pion) propagators (~ (p 2 — m 2 )^ 1 ~ [i?~ 2 ]); 
o fermion propagators (~ (j> — M^)^ 1 ~ [i? -1 ]); 
o derivatives in vertices {p^ ~ [E]); 

o cluster ^-functions ~ [-E~ D ] (C — 1 of them; one has already been removed) ; 
o ^-functions ~ for noninteracting nucleons (see above and App. A); 

an optional factor [in brackets below] (because it depends only on A, and 
not the process) that enforces the overall normalization of the wave function 
in momentum space (see Appendix A). 
In that order we have 

v = DL-2I B -I F + Y^ d t , - D{C — 1) — (D — 1)X N + [(D- 1)(A- 1)] . (6.10) 

i 

There are 5 basic internal "variables" in Eqs. (6.7 - 6.9) (together with 
C, En,V) that determine v(Ib, If, L,bi, fi) and only those 3 relations among 
them; Ib and Ip are always eliminated. We must keep 2 of the remaining ones 
and choose L and The other choices produce equivalent, though less useful, 
formulae (see Ref. [jl2| for a discussion of options). Eliminating Ib, If and J2i 
using Eqs. (6.7 - 6.9) we find 

v = l + (D-2)(n c + L) + J2A-((D-l)(A-l)) + [(D-l)(A-l)] 

i 

= 1 + (D - 2)L + A. (6.11) 

The constant term in square brackets cancels an identical term that arises from 
algebraic manipulation. Note also that the last two terms in Eq. (6.10) in effect 
add (D — 1) (i.e., 3) powers (from phase-space factors) for each independent 
(i.e., non-center-of-mass) nuclear coordinate involved in an interaction. These 
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powers (as we will see below in an explicit example) in effect convert from 
momentum to configuration space. Finally, for D = 4 we have a very simple 
and elegant power-counting result: 

v=l + 2L + A. (6.12) 

Except for irrelevant constant terms (depending on A) , this agrees with Wein- 
berg and van Kolckjll], |i"2f . 

Recall that L and A are separately positive semi-definite. The case L = A = 
[y = 1) cannot occur because this would imply a single nucleon interacting 
without a loop; only the nucleon kinetic energy has this form and it corresponds 
to A = 1, as we found in Eq. (3.7). Thus, the minimum values of v are 2, 
corresponding to the kinetic energy, and 3, corresponding to NN forces (n c — 
1, L = 0, A = 0). See also Appendix B for another example. 

We summarize this section by noting that: 

• Powers, v, of a generic energy or momentum, Q, can be counted in Feyn- 
man diagrams by following the flow of momentum through each vertex. 

• The nuclear case requires consideration of all nucleons, since all nucleons 
eventually interact in a bound state. 

• The final formula for Q v is exceptionally simple, and shows that v always 
increases as more complicated mechanisms are considered. 

7 Electromagnetic Interactions 

We can extend these results to EM interactions within strongly-interacting 
systems by including the "photon" as an extra boson ]T3|. This produces only 
one significant change. For EM interactions: Af M > —1, as noted by RhojLO). 
This is compensated by a factor of e, the fundamental charge, which will reduce 
the size of amplitudes. For specificity we will illustrate the case of a single 
photon (virtual or otherwise) and this will allow us to discuss nuclear EM 
currents (either impulse-approximation or meson-exchange), J EM = (p, J). We 
refer the reader to Refs. Jl(], |2^] for more complex cases. 

Separating out the (assumed) single EM vertex from Eq. (6.12), we have|l(| 

v E m = 1 + 2L + A st + A E m, (7.1) 

where A$t is the sum over strong vertices (> 0) and Aem refers to the single 
EM vertex (> — 1) and will differ for p and J. 

Various important building blocks are shown in Fig. (6). The impulse- 
approximation current, Ji mp , is shown in Fig. (6a). The relativistic form of 
the four-current is N^^N, and involves no derivatives. We know, however, that 
when "odd" 7-matrices are rendered to nonrelativistic form, they are of lead- 
ing order (p/M^v), while "even" ones are 0(1). Thus, Aem = — 1 for Pimp 
(/ = 2,d = 0), while Aem = for 3 lmp (f = 2, d = 1). The pion EM current 
in Fig. (6b) has a single derivative [pi] and Aem = — lj while the dominant 
seagull term^TJ, [3l[] shown in (6c) has Aem = —1 for Jsg an d Aem = for 
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(a) (b) (c) (d) 

Figure 6. Lagrangian building blocks for EM interactions of pions and nucleons. 
Solid lines are nucleons, dashed lines are pions, and wavy lines depict photons. 

Psg- The short-range two-body three-current in (6d) has / = 4, d > 1 and con- 
sequently Aem > 1- This suppresses the short-range contributions compared 
to pion-range and impulse-approximation currents. We note that the Z\-isobar 
current (assuming that we choose not to include the A in our active Hilbert 
space) and the (pttj) current would be (higher-order in 1/A) contributions |Tc| 
of the form illustrated in Fig. (6c), while the (uaj) current fllOf would have the 
form in Fig. (6d). 

We summarize this section by noting that: 

• The electromagnetic current in a nucleus can be treated in a fashion 
similar to the energy, and this is most conveniently done by separating the EM 
vertex from the strong ones. 

• Short-range MEC are suppressed compared to pion-range MEC. 



8 Final Counting of Powers 

Before completing the power counting, it is useful to interpret our results using 
a simple, familiar example. The one-pion-exchange amplitude in Fig. (la) has 
the form 

vAi) = - 9 it tl - t2 Z q ?' q (B-D 
f4 q + m i 

in momentum space, where q is the three-momentum transferred between nu- 
cleons "1" and "2". In terms of power counting, K-(q) ~ Q°/f%. As we argued 
in the previous section and in Appendix A, this should be multiplied by the 
phase-space factor , 2 ? 9 and the Fourier transform to configuration space com- 
pleted. The phase-space factor behaves as Q 3 , in agreement with Eq. (6.12) (for 
n c = 1, L = 0, A = 0). Moreover, performing the momentum integrals produces 
a Yukawa function multiplied by the familiar factor of 1/Att (see Eq. (3.6)). The 
latter converts a single / w in Eq. (8.1) to a A (47r/ w ~ A). Thus (r) ~ -^j- as 
we found in Eq. (3.8). We expect that more complicated diagrams with more 
propagators (n c > l)[fll) will generate a factor of l/(47r)™ c . 
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The interpretation of the short-range force Vsr in Fig. (3d) is somewhat 
different. We again have a factor of (1//^) from Eq. (5.4) and a phase-space 
factor that leads to a (5-function. We write the latter in the form 



S 3 (r) 



1 5{r) 
47r r 2 



(8.2) 



The factor of l/4n is the same as before, and S(r)/r 2 ~ Q 3 , where Q ~ m v 

is the inverse correlation length that sets the scale of the correlation function. 

Thus, each coordinate-space ^-function counts as ~ Q 3 /47r]4]]]. This is actually 

reduced somewhat because of the repulsive nature of the short-range correla- 
te 3 

tions. Nevertheless, combining everything we see that Vsr also counts as 
as predicted by Eq. (6.12) with L = A = and n c = 1. 




Figure 7. Percentages of accrual of kinetic energy (solid line) , potential energy (short 
dashed line), and probability (long dashed line) within an interparticle separation, x, 
for any pair of nucleons in the triton. 



These concepts are illustrated nicely in Fig. (7), which shows the accrual of 
potential and kinetic energy in the triton. A single pair of nucleons separated by 
a distance x 12 is selected, and the expectation value of 0{x— X\-i)0 is calculated, 
where O = T, V, or 1. These values are divided by the value for x = oo, and the 
percentage accrual is plotted. One sees that the major contribution is between 
1 and 2 fm. Moreover, the short-range part of the potential energy is rather 
modest, starting out repulsive and then yielding to the attractive OPEP. One 
sees in these plots that Q ~ ~ [1.4 fm] -1 is a reasonable value. 

We now put together all of the factors for power counting the energy in 
D = 4 dimensions: Q, f^,A. The Lagrangian scale factors are given in Eq. (5.4). 
We also expect a factor of (47r) _2i from loops (see Appendix B for a discussion 
of this and a possible counterexample) and (47r) _,lc from configuration-space 
propagators (see above). This gives our final result for (E), where E is any 
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irreducible contribution to the energy, in terms of the different scales: 

■)l+2(n c +L) + A pV-j-b 



Q\ nc (Q\ 2L+,lc+A 



(E) 



= Q {fJ [V ■ (8 - 3) 

using f + b= J2i( b i + fi) = 2 ^ + 2 K + L) obtained from Eqs. (6.5 - 6.9), 
and A ~ 47r/ w . When counting for other observables (viz., the T- matrix), the 
number of external bosons explicitly enters the equations pi, and this changes 
the factors that set the energy scales. Our formula reproduces all of the previous 
results obtained using less sophisticated techniques. Specific applications are 
relegated to the next section. 

We summarize this section by noting that: 

• Various nuclear energies behave as Q v / A u ~ n " f%° , with more compli- 
cated mechanisms having larger values of v and n c . 

• This formula includes phase-space factors required for conversion to con- 
figuration space, and incorporates momentum sharing between the nucleons. 

9 Results and Discussion 

The most important aspect of this work concerns the relative sizes of iV-nucleon 
forces. The leading-order force connecting N nucleons will have L = 0, A = 0, 
and n c = N — 1, corresponding to the simplest possible calculation (all others 
with L > 0, A > will be smaller). This produces 

Q2N-1 qN 

(Vnnf) - {UA)N -i - Jn^i . (9- 1 ) 

where we have used Q ~ f v in order to make the final estimate. We can also 
use the additional results of Ref. |l6| , who found that 

(V NN ) - 20MeV/pair, 

(V 3NF ) ~ 1 MeV/triplet , 

(Vwf) < TMeV/quartet. (9.2) 

This geometric decrease of the net contribution of many-body forces is consis- 
tent with Eq. (9.1) and is one of the most important results of Ref. pq] , since it 
confirms the role of chiral symmetry in suppressing many-body forces in nuclei. 

Resurrecting the old formalism of Ref. pc{] (applied to pions interacting 
in a nucleus) also provides some insight into the structure of Eq. (8.3). That 
work was predicated upon performing a nonrelativistic expansion of operators 
using a Foldy-Wouthuysen procedure, constructing nuclear operators using the 
superposition principle, and then performing time-dependent perturbation the- 
ory. The nucleons in such a formalism propagate only forward in time, although 
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mesons go forward or backward. Because the operators refer to the entire nu- 
cleus (e.g., J — y~]f—x Ji for some impulse-approximation vertex, Jj), so do 
the propagators. One simply constructs (E — ~ [E -1 ] in terms of the 



nuclear Hamiltonian, H (details can be found in Refs. |L9|, 20 ). Thus, power 
counting should be exactly the same as we have already derived, although it 
will be necessary to redefine the variables. This configuration- space formalism 
automatically incorporates phase-space factors. 

Typical diagrams are shown in Figure (8), with diagrams from Ref. E(J 
on the left expanded as a set of diagrams on the right, where noninteracting 
nucleons are suppressed for simplicity. Although the "nucleus" diagrams sub- 
sume many distinct mechanisms when expanded into "nucleon" diagrams, all 
of the latter share the same topology specified by the former, and this we wish 
to explore. The cross represents the short-range interaction, while the double 
lines represent a nucleus propagator, or a nucleus wave function, and the large 
dots represent the pion-nucleus vertex, J. Not all diagrams are shown. 

Figure (8a) for a nucleus is equivalent to Fig. (8c). Figure (8b) can be 
expanded into Figs. (8d) and (8e). Because the vertices in the "nucleus" for- 
malism contain all nucleons, expanding second-order perturbation theory leads 
to both pion exchanges (8d) and loops (8e) (only one of the loops is shown, 
and this is calculated in Appendix B). Moreover, we see that all possible or- 
derings are included, and that the "loop-like" appearance of the nucleus dia- 
grams results from the forward-propagating nucleus. We count these diagrams 
as v = 1 + 2L + A, where L is the number of "nucleus loops" . If we wish to 
power-count for these "nucleus" diagrams, we must recall that increasing the 
number of connected nucleons by 1 increases v by 2. Thus, we must count the 
short-range interaction as A = 2. This is a simple rearrangement of our original 
form, moving part of n c into A. Note that if one end of a propagator in a nu- 
cleon loop (e.g., in Fig. (8e)) is detached and reattached to another nucleon (as 
in Fig. (8d)), we lower L by 1 and increase n c by 1, keeping L = n c + L fixed, 
and this is why powers of Q depend only on the combination (n c +L) = L. Both 
types count the same and are subsumed in nucleus diagrams (e.g., Fig (8b)). 
Because different integrals generate different factors of (l/4ir), the counting of 
A factors (which leads to Eq. (8.3); see Appendix B) can be different. In this 
example, however, both mechanisms behave as Q 3 /Af v (in leading order). 

Figure (8f) subsumes the graphs of Figs. (8g) - (8i), all of which have v = 5, 
but differing values of L, n c , C, etc. Note that this set of diagrams includes 
three-nucleon forces, vertex corrections to an AW force (only one of which is 
illustrated in Fig. (8h)) and a "recoil" graph (with C = 2) in Fig. (8i) (which 
is a special problem treated later). The graphs of Fig. (8j) comprise those of 
Figs. (8k) - (8n) plus several others that are disconnected (C = 2), and all have 
v = 5. 

Our final example treats the "infrared" singularities discussed in Section 
3. A typical case is Fig. (8o), which subsumes Figs. (8p) - (8s). Nominally 
the graph has L — 2 and thus v = 5. However, if any of the previous 
graphs are sliced in two by a horizontal line, the pion lines are cut, imply- 
ing a minimum energy at that time of Q ~ E v in those propagators. This also 
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Figure 8. Contributions to the nuclear energy from "nucleus" diagrams on the left, 
broken down into the usual "nucleon" diagrams on the right. Single solid lines are 
nucleons, double solid lines represent a nucleus, dashed lines are pions, while a cross 
depicts a short-range interaction in the nucleus approach. 
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holds for Fig. (80) as long as the horizontal line intersects a pion propaga- 
tor. However, a very different result holds when only the nucleus propagator 
(~ E \ H ~ 1/{Q 2 /Mm) ~ ^=pf, rather than 1/Q) is intersected. In that case 
the energy in the propagator is much smaller. These diagrams are therefore 
enhanced, as we discussed earlier, when the Schrodinger equation is solved. For 
this reason, one chooses to define the Schrodinger Hamiltonian in terms of irre- 
ducible operators, which are then incorporated into the Schrodinger equation. 
The "reducible" diagrams not only count differently (we must add a factor of 
A/Q for each infrared propagator, including those in loops, and adjust factors 
of 4-7T [see App. C]), but they also contribute to the "recoil graph" problem. 

That problem arises because graphs of the type shown in Fig. (lb) depend 
on the energy of the nucleus (Eq) and this is contained in nuclear propagators: 
(E n + H — Eo)^ 1 . Because H and Eq ~ Q 2 / 'A and E v ~ Q, it is conventional 
to expand the propagator in powers of (Eq — H)/E 1T ~ ^. Although (Eq — 
H) acting on an initial or final wave function vanishes, it can also kill the 
propagator in the middle of Fig. (80), leading to an irreducible operator. Thus 
the special feature of reducible diagrams leads to a "freedom" in the form of 
reducible operators that we choose to incorporate into our theoretical structure, 
depending on where we place parts of various operators. Experience has shown 
that disconnected, but overlapping graphs, of the type shown in Fig. (8i) can 
be removed (if desired) using a rearrangement, as can the leading order of 
extended overlapping graphs of the type shown in (8g) and (81). This rather 
technical subject can be reviewed in Refs. jl9|, ^0|, ^5j . The practitioner should 
beware of any graph where a nucleon propagator is not required by kinematics 
to carry an energy ~ E v (as in (8g) or (81)) or graphs that are disconnected 
(as in (8i)). 

As an example of how rearrangement affects the power counting, it has 
been shown that the nominal value of v — 5 (obtained from Eq. (8.3)) can 
be changed to v = 6 (as we already saw in Eq. (3.10)) for the 3NF type in 
Fig. (81), and this also holds for the process in Fig. (8g). The short-range 3NF 
resulting from the interaction of 3 nucleons (Fig. (3e)) has n c — 2, A = 1, and 
hence v — 6. Thus, the leading 3NF can be manipulated into v = 6 (rather 
than 5) by a suitable definition of the nuclear Hamiltonian. Note that there is a 
factor of 1/A associated with the additional factor of Q. These different choices 
are neither right nor wrong; they are a theorist's choice. One must simply be 
consistent. 

Our final examples treat a few of the meson-exchange currents in EM inter- 
actions. Using the rules we developed earlier, the power counting for the four 
(n c = 1) graphs in Fig. (9) for J mec gives 2, 2, 4, 4, respectively, assuming 
that the Zi-isobar is treated as a heavy particle (showing that Aem > 1 for 
Fig. (9d) is left as an exercise [^0|). One can show[[i"o[ pp| that graphs (a) and 
(b) contain a factor of (Att/ 2 )^ 1 , [(c) and (d) have an additional factor of A~ 2 ], 
so that these leading-order MEC behave as Q 2 / f^A ~ Q/A, and should be 
comparable to the impulse-approximation result ~ p/M^ ~ Q/A. This pre- 
cisely conforms to the old "rules of scale" . The isobar and heavy-meson MEC 
are suppressed by an additional factor of (Q/A) 2 . Moreover, this counting is 
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Figure 9. Meson-exchange currents of various types. Solid lines depict nucleons, 
dashed lines show pions, while wavy lines illustrate the EM interaction. 



valid on both sides of the current-continuity equation 



v-j(x) = -i[jr,p(x)] 



(9.3) 



using Eq. (9.1). 

The nonrelativistic charge operator is 0(1), and impulse-approximation rel- 
ativistic corrections are 0(Q 2 / A 2 ). The pion-exchange currents (see Appendix 
A of Ref. g§) are 0(Q 3 /4Trf%M N ~ Q 3 /UA 2 ~ Q 2 /A 2 ), consistent with 
Aem = for the seagull charge operator. Note that Q 2 / A 2 is really the same 
as (v/c) 2 , and this also conforms to the old rules of scale [ pl| p0| . 

Additional special cases are worked out in the literature |l2|. p7| . |30| , |3l| , 
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We summarize this section by noting that: 



• iV-nucleon forces scale at least as fast as (Q/A) 1 ^ -1 , implying that two- 
nucleon forces are stronger than three-nucleon forces are stronger than four- 
nucleon forces ... . 

• Results of recent few-nucleon calculations are consistent with this result, 
which makes nuclear physics tractable. 

• The topology of "nucleus" (as opposed to "nucleon") diagrams accounts 
for important aspects of the latter. 

• Infrared singularities (reducible diagrams) enhance the Schrodinger per- 
turbation series. 

• Different treatments of "recoil-graphs" (resulting from IR singularities) 
can lead to changes in power-counting rules (always making terms weaker than 
naive power-counting predictions). 

• Power counting for electromagnetic processes is consistent with the 
current-continuity equation. 

• Heavy-meson MEC are suppressed relative to one-pion-exchange terms. 
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10 Conclusions 



We have developed systematically the power-counting rules of Weinberg, and 
have added additional expository material and numerous examples. Chiral sym- 
metry provides order, and the QCD scale A plays a deterministic role. We 
have shown that the nuclear kinetic and potential energies ( "intrinsic" /pair) 
scale roughly as Q 2 /A, consistent with weakly-bound systems. iV-nucleon forces 
are suppressed as N increases. Increasingly complex contributions to the force 
progressively weaken. Short-range meson-exchange currents are weaker than 
pion-range currents, which are comparable to impulse-approximation currents. 
Large strong-interaction coupling constants of heavy-mesons to nucleons result 
from the mismatch of the scales /„• and A. 

All of these results have been stated before in one form or another using a 
variety of arguments or empirical observations, but their totality rests on power 
counting. In the words of S. Weinberg Ej), 

"The chiral Lagrangian approach turns out to jus- 
tify assumptions (such as assuming the dominance 
of two-body interactions) that have been used for 
many years by nuclear physicists ..." 
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12 Appendix A Momentum Space to Configuration Space 

We have relegated a tedious but instructive part of our derivation to this Ap- 
pendix. We have chosen the option of power counting in configuration space. A 
straightforward power-counting derivation in momentum space does not involve 
the last two terms in Eq. (6.10), which arise from the conversion to configu- 
ration space. These terms are vital, since they reset the baseline (for fixed A) 
against which we determine the importance of various operators. As we will 
see below, it amounts simply to incorporating phase-space factors for each in- 
dependent degree of freedom in a nucleus. In D — 1 space dimensions, there 
are (D — 1)(A — 1) independent internal degrees of freedom, plus D — 1 that 
specify the motion of the nuclear center-of-mass (CM) and don't play any role 
in our discussion. 

We wish to calculate (<^/|0|$i) in both momentum and configuration spaces 
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(for D = 4). We write 

where our notation "{}" emphasizes that the wave function depends on the 
set of internal coordinates {r^} = (r^, r' 2 , . . .) or momenta {p'j}, and we have 
removed the CM coordinates (J2j r 'j = an d J2j Pj = 0). It is a convenience to 
treat the coordinates r'- as independent and use a 8- function as the constraint. 
As a check, we evaluate the normalization integral: 

(3^) ee J n t {<Pr>) 8 3 ^£ r'JAj $} ({r$}) *i({r$}) 

- / (fP^ 3 (EP*) *<«Pi» = 1 • ^ 

The factor of A in the r' k (5-function is conventional. Thus, in order to obtain 
the expectation value of an operator in configuration space, (O), we need to: 
(1) look at the operator in momentum space and (2) multiply by the requisite 
number of independent phase-space factors as given in Eq. (A2). This adds 
[{D — 1)(A— 1)] momentum factors to v (accounting for the last term in Eq. 
(6.10)) and resets the baseline for each diagram or process. Note that if we eval- 
uate ($/|T|!?i), where T is the kinetic energy, we simply insert ^ - p^- 2 /2Mjv 
between the wave functions in the second form of Eq. (A2), which is an obvious 
result. 

If we take the expectation value of a two-body operator, Vijfaj), then a 
somewhat different form results 




f)(p + q'i, -P + q' 2 , • ■ •) V 12 (p) !Pi(q£, q' 2 , . . .) . (A3) 
In addition to the conventional offset ((D — 1)(A — 1)), the two-body poten- 

,3 

tial inherits a phase-space factor: , 2 ? 9 . This is already included in Eq. (6.11) 
because of momentum sharing. The phase-space factors above serve to kill all 
the (5-functions from non-interacting nuclcons, leaving phase-space factors (in 
effect) only for interacting nuclcons. This accounts for the j^Pf m (^3). 
Thus by including momentum sharing with "non-interacting" nucleons and a 
full set of phase-space factors, we have reset the baseline so that our power 
counting works in the same fashion for any nucleus and any operator. Our final 
results do not depend on A. Naive power counting of the potential in momen- 
tum space produces Q° for OPEP; the phase space factor makes this Q 3 , as we 
derived earlier (n c = 1, A = L = — > v = 3). Three-nucleon operators pick up 
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an additional phase-space factor, and so on. This completes the interpretation 
of various factors in the power counting. 

We summarize this appendix by noting that: 

• Inclusion of Fourier-transform phase-space factors resets the baseline for 
all diagrams. 

• This generates a diagram-independent offset [(D — 1)(A— 1)] that makes 
our final power-counting formula independent of A. 

13 Appendix B 



(a) 



(b) 



(c) 



Figure 10. A single nucleon interacting with pions. Solid lines are nucleons, dashed 
lines are pions, and the cross depicts a mass counter term. 

A wide variety of concepts can be illustrated by working out a simple ex- 
ample in toto. Figure (10a) illustrates the self-energy of a nucleon arising from 
pionic vacuum fluctuations. We perform a nonrelativistic calculation of the 
S-matrix, noting that the nucleon self energy ^2 N is traditionally defined by 
S = -iE N [(2n)W(P f - Pi)\- 



J N 



d 4 k 



k)H 2 



(2vr) 4 {k l-El + l e)(-k Q + E-^f +l e) 



(Bl) 



using the Feynman rules: ^t a cr ■ k for creating a pion with isospin component 
a, momentum k^, and spin er, which corresponds to A = 0. We have included 
(for now) the complete nucleon propagator (energy E = p 2 /2M/v, momentum 
p), and have defined the pion energy by i? 2 = k 2 +m 2 . Using (<r-k) 2 = k 2 , t 2 = 
3/4, and noting that there is a single pole in the lower Aft-plane we obtain 



8 \.U 



d 3 k 



(2tt) 3 



& 



'■IT (yE v 



(p + k)2-p- 

2M N 



r 



(B2) 



To leading order in 1/Mjy, we can ignore the kinetic-energy terms ~ 1/Mjv. 
We also note that the pion energy sets the scale in both of the propagators. 
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Much of the power counting that we did earlier depends explicitly on this fact. 
Thus to leading order in (1/Mn), we obtain the cubically divergent result 

Interpreting this result requires further work. First, we must regularize the 
integral and render it finite. Then we must renormalize. The method of choice is 
dimensional regularization. We convert the 3-dimensional integral to n "space" 
dimensions (J , with n not necessarily an integer), and evaluate it for 
some n where Em is finite, renormalize the result, and analytically continue 
this finite value to n = 3. The angular integrals plus factors of 2tt give[[45| 

<Mn = 2 

(2tt)« r(n/2)(47r)«/2 ' ^ > 

and we see the origin of factors of 47r. The remaining integral isQ] 

dk k n+1 ml TT 



(fc 2 +m2) 2 sin(2f) 



(B5) 



For n < these integrals are finite, and their extension to n = 3 (or any odd 
dimension) is also finite. Combining Eqs. (B4) and (B5) for n = 3 produces 
m^/47r. 

At first sight this is a very strange procedure, and the reader is referred to 
Ref. Q for more expert justification. We note that by adding and subtracting 
m\ in the numerator of the integrand in Eq. (B3) we obtain J ^~ft 9 , which has 
no length or energy scale and therefore has no obvious physics associated with 
it (so we drop it), plus another term. This process can be repeated, producing 

another scaleless integrand (— m\ J p^prp-) plus ^ , the last part of which 

we obtained previously using Eqs. (B4) and (B5). Thus, we finally obtain a 
well-known result |L4j 

Zg\ml 3g 2 A ml m 3 n 

N 327T/2 8U(4ttU) ~ UA ' [ 0) 

Manipulations of the type used here to interpret our results should never be 
performed when they introduce singularities at k = (ours were finite there). 

What is the interpretation of our procedure, and why did we keep a single 
term and argue away the rest? Figure (10b) also occurs as a part of Em- This 
term is simply M N , the "bare" nucleon rest mass. The divergent (but scaleless) 
integrals can be viewed as contributions to M N . Since we did not know M N 
anyway, adding terms to it does not introduce a complication. The same is 
true for Em above, but its properties are special. Our original methodology 
was to divide the physics into soft (long-range) and hard (short-range) parts. 
The scaleless terms are hard (divergent), but Em is soft. It is proportional 
to m\ ~ Q 3 in accord with our counting rules (L = 1, n c = 0, A = 0) and 
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is nonanalytic in m^. That is, it has the form y/m^ = rrv^yfm^. Because 
pion masses must appear in the Lagrangian as powers of m 2 , the square root 
is special. Nonanalytic terms are usually logarithms, but not always. It is a 



common practice^, 34 to separate out all nonanalytic (soft) terms and lump 
all analytic terms (even finite ones) with the so-called "counter" terms (hard 
terms such as Mjy). Nevertheless, because the nonanalytic terms have a special 
form (and are often large), we make the separation: Mm = M N + En- Note 
that Mat ~ A, but En ~ m%l fnA ~ rn^/A and is significantly smaller, as we 
expect on the basis of power counting (L = 1, n c = 0, A = — * v = 3). 

Finally, we note that a single factor of l/4ir arose in Eq. (B6). This is 
somewhat unusual, as normally loop integrals generate l/(47r) 2 or something 
similar. A single (l/47r) arises in nonrelativistic cases corresponding to odd- 
dimensional (for n = 3) integrals, such as Eq. (B3). In each case the nonanalytic 
term is an odd polynomial. In the usual case such as Fig. (10c) one obtains 
l/(47r) 2 and a logarithm. A useful and instructive exercise is to calculate Fig. 
(10a) covariantly. One finds that, indeed, all terms generate an explicit factor 
of l/(47r) 2 . If one expands that result for small /Mm , one finds 

rRei _ aMff + b M N m\ + c mj/M N 3g>4 log(m„/M N ) 

where a, b, and c are dimensionless, and a and b are divergent. The a — c terms 
are analytic functions (polynomials) of m 2 . and can be incorporated directly 
into Mjy, as we discussed above. The two scaleless integrals that we found 
above contribute to a and b, respectively. The logarithmic term is unique in 
the expansion. 

In performing the expansion in Eq. (B7) that produces En (Eq. (B6)), 
a factor of tt/2 is generated, leading to a single residual factor of tt in the 
denominator of Eq. (B6). The logarithmic term has a dimensionless coefficient 
of 3g\/2 ~ 2.4 and a factor of l/(47r) 2 , which was assumed in our derivation 
of Eq. (8.3) that leads to E N — Q 3 /A 2 . If wc force Eq. (B6) into this form a 
very large dimensionless coefficient of —3^3^/2 ~ —7.5 results. Sometimes this 
happens. 

Our nonrelativistic calculation reproduced the leading-order nonanalytic 
part of the covariant calculation (and was much easier) . Although the analytic 
parts of the two calculations are different, they are not required to be the same 
and this cannot affect the final results, since we do not a priori know M^. 

We summarize this appendix with the following observations: 



• Sensible nonrelativistic field-theory calculations are possible. 

• Dimensional regularization is an easy way to make integrals finite. 

• Loops can be rendered into analytic (typically "hard") parts that are 
incorporated into coupling constants and nonanalytic ("soft") parts, which are 
kept separate. 

• Power counting works for loops. 

• The pion mass controls the scale of loop propagators. 
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• Factors of l/(4-7r) 2 (and occasionally (l/47r)) arise from loops. 

• Since A ~ 47r/ 7r and Q ~ /„. , counting powers of Q /.4 is not very different 
from counting powers of (l/Air). 

14 Appendix C: Zero-Range Model 

Another excellent example is the zero-range force pi). We write the Lagrangian 



for two identical nonrelativistic nucleons|47 interacting via a zero-range force 
as 

C = iV+(x,t)(<^ + J^-)JV(x,t) - A(7Vt(x,i)7V(x,t)) 2 . (CI) 

The form of the scattering amplitude generated by this Lagrangian is a series 
of loop diagrams, each order (in A) being a product of loops involving the 
two nucleons. Performing the ko integral (part of the d A k/(2n) A phase-space 
factor in each loop) leads directly to the Schrodinger equation corresponding 
to an energy E = fc 2 /2/i and the nuclear (CM) Hamiltonian:p 2 /2/i + 2A<5 3 (r), 
where r is the separation of the nucleons and /i is the reduced mass. We have 
combined the kinetic energies of the two nucleons so that Mn is twice /i; note 
the combinatorial factor (of two) in front of the potential. That equation can 
be written in the formQ (note that G (r) = -§f ^): 

!P(r) = !P (r) + f d 3 r' G (r - r') V(v') V (r') , (G2) 

where V(v') is the (J-function) potential, !?o( r ) is a plane wave, and Go(z) is 
the Green's function: (E — H + ie) _1 . Performing the integral, we obtain 

ff(r) =^o(r) + 2AG (r)^(0), (G3) 

the "zero-range" solution; setting r to produces an algebraic equation 

from which we obtain the T-matrixp7| (unitarity implies Im(T) = 
-2^|T| 2 /47r): 

T= 1-2AG (0)- (C5) 

The Green's function at the origin is (of course) linearly divergent, but we can 
apply the regularization techniques of Appendix B: 

G °<°> = 2 */ » {ce) 

—kfii 
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where again the "rt-dimensional" integral is finite (and is left as an exercise, 
being only slightly different from Eq. (B5)). This produces 

t=y^m, (CI) 

7T 

which corresponds to a scattering length, a = fj,\/ir, an S-matrix, 

S = (1 — ika)/(l + ika), and an effective-range function (the inverse of the 

K-matrix), k cot(<5) = —1/a. 

We know that poles in the T-matrix indicate special states. A single pole 
always exists in this model for 

k 7F * 
/iA 

and 

E = ^~l/X>. (C8) 

This corresponds to a bound state if A > and to a "virtual" state if A < 0. 

This is a most peculiar result, since we started with a potential that is 
repulsive if A > 0, and as A — > the bound state gets deeper! One must remem- 
ber that the original problem (interpreted as a quantum mechanics problem) 
has no solution at all. We have generated a solution by changing the prob- 
lem, redefining Go(0) and making it finite, and the peculiar properties of (C8) 
are a reflection of the original (insoluble) problem. This redefinition is equiva- 
lent to defining how loops (vacuum fluctuations) contribute. Thus, A is really 
the renormalized coupling constant (of arbitrary sign after renormalization), 
and not the "bare" one in Eq. (CI). How one might interpret and treat these 
peculiarities is discussed in Ref. psfl . 

What about A — » ? This won't happen in general, since applying Eq. (5.4) 
leads to the equivalence 

A^, (C9) 

Jit 

with |ca| ^ 1, and using 2fx — Mn ~ A ~ 47r/ 7r , we find 

K = ^4^~^r~Q, (cio) 

(2/x)2ca 2c A 

- f 2 

as well as E ~ . j T . and a ~ 2c\/f n . A bound state with binding energy ~ 2.5 

A 

MeV and a corresponding scattering length ~ 4.3 fm is generated for c\ ~ 1. 
The case c\ ~ — 1 produces a scattering length ~ —4.3 fm. We note that an 
anomalously small ca is very improbable. 

The unitarity of the T-matrix in lowest-order PT relates the first- and 
second-order amplitudes with a single factor of l/(47r). The (second-order PT) 
loop integral must therefore be of the type that generates a single factor, rather 
than the usual l/(47r) 2 . 
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